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ABSTRACT 

We present the first spectroscopic measurements of the [Oiii] 5007A hne in two z ^ 3.1 Lyman- 
alpha emitting galaxies (LAEs) using the new near-infrared instrument LUCIFERl on the 8.4m Large 
Binocular Telescope (LBT). We also describe the optical imaging and spectroscopic observations used 
to identify these Ly-a emitting galaxies. Using the [O ill] line we have measured accurate systemic 
redshifts for these two galaxies, and discovered a velocity offset between the [O iii] and Ly-a lines 
in both, with the Ly-a line peaking 342 and 125 km s""'^ redward of the systemic velocity. These 
velocity offsets imply that there are powerful outflows in high-redshift LAEs. They also ease the 
transmission of Ly-a photons through the interstellar medium and intergalactic medium around the 
galaxies. By measuring these offsets directly, we can refine both Ly-a-based tests for reionization, 
and Ly-a luminosity function measurements where the Ly-a forest affects the blue wing of the line. 
Our work also provides the first direct constraints on the strength of the [O iii] line in high-redshift 
LAEs. We find [O iii] fluxes of 7 and 36 xlO"^'' erg s"^ cm^^ in two z 3.1 LAEs. These lines 
are strong enough to dominate broad-band fiux measurements that include the line (in this case, 
band photometry). Spectral energy distribution fits that do not account for the lines would therefore 
overestimate the 4000A (and/or Balmer) break strength in such galaxies, and hence also the ages and 
stellar masses of such high-z galaxies. 

Subject headings: galaxies: high redshift — intergalactic medium 



1. INTRODUCTION 

The Lyman-a emission line is a highly efficient tool 
for identifying and studying star forming galaxies at 
high redshifts. This line can carry up to 6% of the 
bolometric luminosity of a young stellar population 
([Partridge fc PeeblesI il967t ) and is conveniently placed 
for observations by ground-based optical observatories 
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for 2 < z < 7. However, the transmission of Ly-a emis- 
sion is complicated by its resonant scattering interaction 
with neutral hydrogen, both within the galaxy emitting 
the line and in the surrounding intergalactic medium 
(IGM). The Ly-a line is observed in abou t 25% of z ^ 
3-5 Lyma n-break galaxies (LB Gs) fe.g.. [ Steidel et aTl 
l2000t iDo^Hygelund et ah .2007c iRhoads et al.ll2009D . a 
percentage that may increase with incr easing redshift 
([Shimasaku et al.ll2006l: IStark et al.ll2010t) . This trend is 
supported by work at redshifts less than two, where the 
fraction of galaxies exhi biting Ly-a emission decreases 
at these lower redshifts ([Reddv et al. 2008; Hayes et al] 
120101: iCowie et al.]|2010f ). The Ly-a line is observed to 
have a characteristically asymmetric profile, with a sharp 
cutoff on the blue side and a more extended wing on the 
red side fe.g.. IRhoads et al]l2003D . In Lyman-break se- 
lected galaxies, the peak of the Ly-a line is typically 
redshifted by several hundred kms~^ with respect to in- 
terstellar absorption l ines and/or nebula r emission lines 
([Shaplev et al.l 120031: iSteidel et ai]|2010[ ). whereas this 
measurement has not been made in Ly-a selected galax- 
ies until this paper. 

Besides being a useful tool for studying galaxy proper- 
ties, Ly-a galaxies also offer unique and powerful probes 
of cosmological reionizatio n (e.g., [Rhoads & M alhotra 



2001 1 : iMalhotra fc Rhoads|[200H2006t iKashikawa et al 



20061: IMcCandlissI 120091: IPaval et alJ iMoi) The de 

tailed interpretation of these tests can be substan- 
tially affected by velocity offsets between Ly-a and 
the systemic redshift, because a redshifted line is 
less affected by the damping wing of Ly-a absorption 
from the IGM (ISantosI 120041: IMalhotra fc Rhoadsl [20061: 
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IDiikstra fc WvithdIMol) . 

It is not sufRcient to assume that the velocity off- 
sets seen in LBG samples hold for Ly-a selected sam- 
ples. LAEs are typically less massive than presently 
available Lyman-break selected samp les at s imilar 
redshifts (iVen emans ct al., ,2005u .Gawiser et alj 120061: 
Finkelstein et a l. 2007; Pirzkal et al.ll2007l : iNilsson et alJ 
20071 ). They should have correspondingly lower escape 
speeds, provided that stellar mass correlates broadly with 
dark matter halo mass. Such a trend is discussed by 
[Gawiser ct al. (2007,), where their sample of z ~ 3.1 
LAEs have typical stellar masses of 1 x 10^ Mq and 
median halo masses of 7.9 X lO(lO) Mq. Gawiser et 
al. point out that these values are significantly smaller 
than those values for LBGs at z ~ 3.1, which have 
stellar masses of '-^ 2 x 10^° Mn and halo masses of 
-3 X 10^1 Mq (jShaplev et al.|[200ll : lAdelberger et al.l 
I2005f ). Galactic winds (and indeed many other astro- 
physical outflows) typically have flow speeds near the 
escape speed for the object, and the observed velocity 
offset of a Ly-Q line is roughly double the wind speed 
(jVerhamme et al]|2006[ ). Additionally, the velocity off- 
sets in Lyman-break samples are inversely correlated 
with the Ly-a e mission strength, as characterized by 
equivalent width (jShaplev et al.l[200ll ). and the equiva- 
lent widths of the Ly-a selected samples are much larger 
on average than those of LBG samples. Finally, Ly- 
g selected galaxies are typically small in ph ysical size 
(|Bond et al.l[20Ql[20T0l:[Malhotra et al.l[2010l) . 

We present here the first direct measurements of the 
velocity offset between Ly-a and nebular emission lines 
for typical Ly-a selected galaxies. Our measurements 
are based on a combination of near-infrared spectroscopy 
with the new LUCIFER instrument on the Large Binocu- 
lar Telescope, and optical spectroscopy using Hectospec 
on the MMT. We selected targets for the study from 
a large area narrowband survey conducted with the 
OOPrime camera on the 2.3m Bok telescope of the Stew- 
ard Observatory. 

In section [5] we describe our observations and data 
analysis methods. We present our observational results 
in section [31 and discuss their implications in section |4l 
Where relevant, we ad opt Hp = 70 km s ~^ Mpc^^, 
Um = 0.3, flA = 0.7 (iSpergel et al.l I200I . Also we 
use the following vacuum wavelengths, 1215.67 A for 
Ly-a, 3729.875 A for [On], 4862.683 A for and 
4960.295/5008.240 for [Oiii] from the Atomic Line List 
v2.04E|. All magnitudes quoted are AB magnitudes un- 
less otherwise specified. 

2. OBSERVATIONS AND DATA PROCESSING 

2.1. Narrowband Survey - Observations and Data 
Reduction 

We completed a deep narrowband survey for LAEs at 
z ^ 3.1 using the 90Prime Camera on the 2.3m Bok tele - 
scope at the Steward Observatory (|WiUiams et al.|[200^ . 
The survey was completed in the COSMOS field cen- 
tered at RA 10:00 :28.6 and DEC -1-02:12:21.0 (J2000) 
(|Capak et all [2001 . The KPNO MOSAIC [Oiii] filter, 
centered at 5025 A, with a bandwidth of 55 A, was used 
to select Lyman alpha emission at redshifts z ^ 3.1. 

http:/ /www.pa.uky.cdu/'^pcter/atomic/index.html 



The data was obtained through time allocated by Stew- 
ard Observatory in February 2007 (PI Finkelstein) and 
February and March 2009 (PI McLinden). We have cre- 
ated a 1.96 deg^ image, representing a total integration 
time of 16.67 hours. The complete details of this survey 
and the data reduction process will be highlighted in a 
forthcoming paper. 

2.2. Broadband Data 

We obtained publicly available broadband imaging 
data in CFHT u* and SDSS g+ bands from the 
NASA/IPAC Infrared Science Archivj^ to complement 
our narrowband survey. The g^ im aging data (v2.0) 
comes from Suprime-Cam (jMivazaki et al. 2002.) on the 
Subaru 8.3m telescope. The u* band imaging data (v5.0) 
comes from the MegaPrime/MegaCanQ on the Canada- 
France-Hawaii 3.6m Telescope. The 5cr depth in a 3" 
aperture in each ban d is 26.4 and 27.0 f or the u* and g+ 
bands, respectively (jCapak et al.l 120071 ). These broad- 
band images were registered to our narrowband image 
using the IRAF tasks WCSMAP and GEOTRAN, which 
resamples the broadband images to match the coordinate 
system of the narrowband image (0.45 "/pixel). 

2.3. Candidate Selection from Narrowband and 
Broadband Data 

We used the SExtractor package (jBertin fc Arnout^ 
11996') to perform source detection in the narrowband and 
broadband images. SExtractor was run in dual-image 
mode, first with narrowband image as both the detec- 
tion and measurement image, and a second time with 
the narrowband as the detection image and broadband 
image as the measurement image. We selected LAE can- 
didates based on the strength of their narrowband versus 
broadband excess as well as their colors as outlined in 
Rhoads & Malhotra (2001). 

Namely, LAE candidates must be detected in the nar- 
rowband at the 6(7 level, their flux in the narrowband 
must exceed that in the broadband (g+ band) by at least 
a factor of 2 and their narrowband flux must exceed their 
g^ band flux at the 4(7 level. Finally, candidates must 
have flux in the filter bluward of of the Ly-a line (u* 
band) consistent with expected Ly-a forest absorption 
blueward of the Ly-a line and consistent with a u* - g^ 
color > 2. Selection criteria are shown below in equations 
1 - 4. 

fnb/Sfnb>G (1) 

W/s>2 (2) 

fnt-f,>HSf^, + Sf^f (3) 

fu < 10-^/^5 + (3 X SU) (4) 

Here inb is the narrowband flux, fg is the g band flux, /„ 
is the u band flux, Sfnb is the flux error in the narrow- 
band, Sfg is the flux error in the g band, and Sfu is the 

http;/ /irsa.ipac.caltech.cdu/data/COSMOS/datascts.html 
Based on observations obtained with MegaPrime/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the Na- 
tional Research Council (NRC) of Canada, the Institut National 
des Science de I'Univers of the Centre National de la Recherche 
Scientifique (CNRS) of France, and the University of Hawaii. 



3 



flux error in the u band. The second criterion requires 
that objects have Ly-a equivalent widths > 57.5 A. The 
decision to require a 6cr detection in the narrowband and 
to require a 3a non-detection blueward of the Lyman- 
break indicates that these are stringent selection crite- 
ria, meant to exclude false detections and low redshift 
interlopers. We used isophotal magnitudes (MAGJSO) 
from SExtractor to measure the magnitudes and fluxes of 
each object. Isophotal magnitudes were chosen because 
they have been found to produce the most accurate colors 
when SExtractor is run in dual-image mode (Holwerdal 
I200I) . Isophotal magnitudes are not measured within a 
fixed aperture for each object, but rather determines the 
magnitude from the number of counts in pixels above the 
user defined threshold and hence each object has a unique 
'aperture' in which its flux is measured. For instance, the 
two LAEs with detected [O iii] emission have extracted 
narrowband isophotal areas of 20.86 arcsecond^ and 8.71 
arcsecond^ (later referred to as objects LAE40844 and 
LAE27878, respectively). Similar selection criteria used 
at z=4.5 have typically yielded a spectroscopic success 
rate of 80% (Dawson et al. 2004, 2007, Rhoads et al. 
2003, 2005, Wang et al. 2009). 

2.4. Optical Spectroscopy - Observations and Data 

Reduction 

We obtained optical spectra of LAE candidates in 
January, February and April 2009, using the Hectospec 
multi- fiber spectrograph at the 6.5m MMT Observatory 
(a joint facility of the Smithsonian Astrophysical Obser- 
vatory and the University of Arizona). Hectospec has a 1 
deg^ field of view and spectral coverage from 3650 - 9200 
A. The resolution of the instrument is 6 A. Optical 
spectra are crucial for confirming that candidates are in 
fact LAEs at the correct redshift and not lower redshift 
interlopers and for determining the exact wavelength for 
the Ly-a line. 

We reduced the Hectospec data and extracted ID spec- 
tra using the External SPECROAEF^ pipeline developed 
by Juan Cabanela. The External SPECROAD pipeline 
applies bias, dark and fiat field corrections as well as 
wavelength calibration using He-Ne-Ar arc lamps. Typ- 
ical residuals from the wavelength calibration are 0.15 

A. 

The optical spectra of the three objects chosen for 
near-infrared followup are shown in Figure 1. These 
optical spectra confirm that these objects are in fact 
LAEs. The spectra show strong Ly-a emission at the 
expected wavelength and the line displays the charac- 
teristic asymmetry e xpected for this line when emitted 
from a high-z source ([Rhoads et al. 2003; Dawson " et al.l 
[200I IKashikawa et al.ll20'06ir See section 3.1 for further 
discussion of this asymmetry. Finally the spectra were 
checked for the presence of any other optical lines. No ad- 
ditional emission lines were observed at the wavelengths 
where they might be expected for foreground [0 11] or 
[O III] emission line objects. 

2.5. Near-Infrared Spectroscopy - Observations and 

Data Reduction 




Fig. 1. — Hectospec optical spectra - used for confirmation of 
the objects as Ly-o emitting galaxies at z ~ 3.1. LAE40844 has 
an additional feature, a 'blue bump', at A ~ 4990 A, which is 
discussed in section 3.4. The spike seen near 4982 A is noise and 
not an additional feature. 

Three of our brightest confirmed LAEs after Hectospec 
observations were observed in the near-infrared (NIR) us- 
ing the new near-infrared instrument LUCIFER (LET 
NIR Spectrograph Utility with Camera and Integral- 
Field Unit for Extrag alactic Research) on the 8.4m 
LET (Seifert et al. 2002i; lAgeorges et al.l[20Tol ) The Ly-a 
line flux of these objects chosen for NIR followup, de- 
rived from their narrowband and broadband magnitudes, 
ranges from 0.94 - 3.6x10-1^ erg s^^ cm'^. LUCIFERl 
is the first of two planned NIR instruments for the two 
8.4m mirrors of the LET. LUCIFERl currently operates 
on one mirror of the LET and is capable of spectroscopy 
and imaging in the wavelength range 0.85/Ltm - 2.5/im. 
Our observations were performed in queue mode during 
LUCIFER'S Science Demonstration Time in December 
2009 and continued during the Arizona Queue in Jan- 
uary and February 2010. 

We used the longslit mode of LUCIFER with a 1" 
slit utilizing the H-t-K grating with 200 lines/mm and 
the N1.8 camera. The image scale of the N1.8 camera 
is 0.25"/pixel. We obtained 10 two- minute integrations 
for our brightest LAE (henceforth LAE40844) . Our sec- 
ond object (henceforth LAE27878) was observed over 20 
four-minute integrations. Our final object (henceforth 
LAE14310) was observed over 25 four-minute integra- 
tions. The central wavelength in this setup is 1.93 /xm, 
and the spectral coverage spans essentially the full H and 
K band windows. The spectral resolving power with the 
4 pixel slit ranges from 940 near 1.6 fim to 1286 near 
2.2/im or a resolution of ~ 4.3 A/pixel. 

We utihzed the DOSLIT routine in IRAF (jValdesI 
Il993f ) to reduce the 2D spectra. To simplify reduction, 
our observations were designed so that a bright (R 12 
- 18) continuum source also shared the slit with each 
LAE. This allows for a trace to be created using the 
bright continuum object. The trace was then shifted 
along the spatial axis to extract the much fainter LAE, 
whose continuum emission is undetectably faint in indi- 
vidual exposures. We performed fiat fielding and dark 
correction before aperture extraction. An aperture of 6 
pixels was used for extraction. Wavelength calibration, 
also performed as part of the DOSLIT task, was done us- 
ing an argon lamp spectrum observed in the same setup 
as our science observations. After reduction, we aver- 

http://iparrizar. mnstate.edu/ -juan/research/ESPECROAD/indexg^ individual frames USing the IRAF task SCOMBINE 
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to produce a single averaged spectrum for each object. 
Average RMS uncertainties from wavelength calibration 
for the two spectra with [Oiii] detection were 0.64 A 
and 0.48 A for LAE40844 and LAE27878, respectively. 
Residual bright night sky lines, a problem when extract- 
ing faint sources, were interpolated over in each night's 
averaged spectrum using the SKYINTERP task found 
in the WMKONSPEC package designed for Keck NIR- 
SPEC reductiorE3- 

Flux calibration proceeded using the bright contin- 
uum sources that shared the slit with our LAEs as de- 
scribed in the paragraph above. Henceforth these con- 
tinuum objects will be called calibration stars. This pro- 
cess corrected for telluric absorption and transformed our 
flux to Fa units. LAE40844 was calibrated using SDSS 
J100126.08+021902.2 and LAE27878 was calibrated us- 
ing SDSS J100025. 10-1-022552.0. We flux calibrated each 
night's calibration star spectrum usin g an appropriate 
Pickles model spectrum (|Pickleslll998l ). scaled in flux to 
match the apparent V magnitude of the observed cali- 
bration star. The appropriate Pickles model was chosen 
based on the spectral type of the calibration star and 
spectral type was determined from SDSS u-g and g-r col- 
ors of the calibration stars as outlined in iFukugita et al.l 
jloT^. The SDSS u, g and r magnitudes come from 
SDSS DR7. The V magnitude of the observed of the cal- 
ibration star was determined from its SDSS colors and 
the Lupton (2005) color transformation from SDSS g-r 
color to V magnitudeEEI. We then created a sensitivity 
curve by dividing the scaled down Pickles model by the 
calibration star's stellar spectrum in counts. We then 
multiplied each night's LAE spectrum (in counts) by that 
night's sensitivity curve to produce a final flux-calibrated 
LAE spectrum. 

2.6. Cross check of photometric redshift 

We cross checked the coordinates of each LAE with the 
sources in th e COSMOS Photometric Redshift C atalog 
Version 1.5 (|Ilbert. Capak fc Salvato et al.l 120091 ). We 
found a unique match for each object, within 1" in all 
cases. The photometric redshift of LAE40844 is Zphot = 
3.094, with a 68% confidence range of 3.08 < Zphot < 
3.11. The photometric redshift of LAE27898 is 3.086, 
with a 68% confidence range of 3.02 < Zphot < 3.12. 
Finally, the photometric redshift of LAE14310 is 3.035, 
with a 68% confidence range of 2.98 < Zphot < 3.11. 

2.7. Cross check with Chandra COSMOS X-ray Sources 

We also compared the locations of our LAEs with 
the C handra COSMOS Survey Point Source Catalog 
(| Elvis et al. 2009) to exclude contamination from AGNs. 
The Chandra COSMOS Survey Point Source Catalog 
contains 1761 X-ray sources in the full 0.5 - lOkeV band 
with a limiting depth of 5.7xl0~^^ erg cm~^. The 
survey covers the central ^ 0.9 deg^ of the COSMOS 
field. We find no X-ray sources matching the coordi- 
nates of any of our LAEs within 12.8", which is much 
larger than the combined positional uncertainties of the 
narrowband and X-ray catalogs. This gives upper limits 
fx/ fhy-a ^ 1-6 - 6.0 for the three sources with LUCIFER 

^'^ http:/ /www2. keck. hawaii.odu/inst/nirspec/wmkonspcc. html 
http:/ /www.sdss.org/dr6/algorithms/sdssUBVRITransform.h 



spectra — below the typical ratio fx/ fhy-a ~ 8 for type I 
quasars, and overlapping the range (fx/ fhy-a ^ 3 - 4 ) 
observed for type II quasars (e.g., iZheng et al.l l20ldl) . 
Thus the X-ray observations suggest that the Ly-a in 
these objects is indeed powered by star formation rather 
than AGN activity, though the present X-ray data are 
not deep enough to prove this case by themselves. Also 
we note that the modest [O ill] velocity widths of ^ 200 
- 300 kms~^ seen in our two LAEs are much lower than 
the typical velocity widths of around 1000 kms"-'^ ex- 
pected for Type 1 AGN. 

3. RESULTS 

3.1. Gaussian Fits to the [Olll] and Ly-a Lines 

To determine the central wavelength and line flux of 
each emission line we fit a Gaussian plus constant to 
each emission line. In the case of the [O ill] line, we fit a 
symmetric Gaussian to the line using the MPFITEXPR 
IDL routine, which is part of the MPFIT packageEB For 
the Ly-a line we fit an asymmetric Gaussian by modi- 
fying the ARM.ASYMGAUSSFIT IDL routineE], which 
also utilizes the MPFITEXPR routine. In its unmodified 
form, ARM_ASYMGAUSSFIT basically fits the left and 
right sides of the central wavelength with different Gaus- 
sians and then requires that in the final fit the left and 
right curves must have the same center and same ampli- 
tude where they meet, meaning there are eight parame- 
ters, four for each side of the curve (amplitude, center, 
sigma, constant) but only six of these are free parame- 
ters. This allows for a single curve to be fit, but the curve 
can have different sigma values for the right and left sides 
of curve, making it ideal for fitting a Ly-a line with a 
truncated blue side and extended red wing. We modified 
ARM.ASYMGAUSSFIT by fixing the constant on the 
left side of the Ly-a emission line to a pre-determined 
constant measured as the average continuum level from 
4000 - 5000 A. This reduces the number of free parame- 
ters from 6 to 5 when fitting the Ly-a line. The constant 
on the right side of the Ly-a line is allowed to vary, as 
one can expect a slightly higher continuum level redward 
of rest-frame Ly-a. 

We quantified the asymmetry of the Ly-a peaks by 
defining the ratio of the redside best-fit sigma to the blue 
side best-fit sigma, or asymmetry = cfred/ <^biue- Using 
this definition, any line with an asymmetry measure > 
1.0 is considered asymmetric. From this definition of we 
find asymmetry measurements of 1.1 ± 0.1 , 2.1 ± 0.2 and 
1.0 ± 0.1 for LAE14310, LAE40844 and LAE27878, re- 
spectively, meaning the Ly-a line in LAE40844 is highly 
asymmetric, whereas the Ly-a lines in LAE14310 and 
LAE27878 appear to be symmetric within the errors. For 
comparison with other asymmetry measurements in the 
literature we also ca lculated asymmetry using a^ and a/ 
(jRhoads et al.l[2003[ ) from the best fit asymmetric Gaus- 
sians. aA is 1.2, 2.2 and 1.2 and a/ is 1.2, 1.8 and 1.1 for 
LAE14310, LAE40844 and LAE27878, respectively. 

We defined the redshift of the emission line using the 
central wavelengths determined from these fits (from z 
= (Aobs/Ae„i) —1 where Aem is the rest-frame vacuum 

developed Craig Markwardt 

|http: / /www. physics. wise. edu/~craigm/idl/idl.html| 

^' developed by Andrew Marble 

http: / /bubble. as. arizona.edu/idl/arm/ 
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wavelength and Xobs is the central wavelength of the best 
fit). Line flux for the [O ill] line was determined from the 
area under the best fit symmetric Gaussian. Line flux for 
the Ly-a line was determined from the narrowband line 
flux. The area under the best flt asymmetric Gaussian 
was scaled to match this flux. The passband of the fil- 
ter transmission curve was taken into account to assign 
the appropriate amount of flux to the main Ly-a peak in 
LAE40844 and this object's secondary 'blue bump' dis- 
cussed in more detail in section 3.4. Errors on the area 
were determined directly from the PERROR output from 
the MPFITEXPR routine, which returns the one-sigma 
error on fitted parameters. PERROR output values were 
also used to quantify the error on the best fit central 
wavelength, but an additional error term was included 
here to account for wavelength calibration errors from 
the Hectospec and LUCIFER spectra. Errors on calcu- 
lated values for redshift and velocity offsets between the 
Ly-a line and the [O ill] line were derived from best fit 
central wavelength and its associated error as described 
directly above. Best fits are calculated from unsmoothed 
spectra, while the spectra in figures 2 and 3 are plotted 
after 3-pixel boxcar smoothing. 

3.2. [Oiii] Detection with LUCIFER 

We detect the [Oiii] 5008.240/4960.295 A doublet in 
two of the three LAEs, LAE40844 and LAE27878. For 
the stronger [O ill] line (rest frame vacuum wavelength 
of 5008.240 A), we measure a line flux of 35.48 ± 1.15 
X 10-" erg s-i cm-2 in LAE40844 and 6.96 ± 0.33 x 
10"" erg s-i cm-2 in LAE27878. 

The second strongest [O ill] line (rest frame vacuum 
wavelength of 4960.295 A) was also found in the same two 
LAEs. The line fluxes measured for this line from best 
fit Gaussians were 14.82 ± 2.24 x lO^^'^erg s"^ cm"^ 
and 1.47 ± 0.37 x lO'^erg cm^^ for LAE40844, 
LAE27878, respectively. The ratio of this secondary 
[O III] line to the stronger [O ill] line is within 2a of the 
theoretical value (1/3) in both galaxies. This provides a 
check of the data analysis and increases confidence that 
this is the 4960.295 A line. Table 1 summarizes the [O III] 
and Ly-a line fiuxes for each LAE, along with relevant 
broadband and narrowband characteristics. 

LAE14310 showed no detectable [Oiii] emission. This 
could be explained if LAE14310 was in fact a lower red- 
shift interloper, but a visual inspection of the optical 
spectrum has ruled out the object as an [O ii] emitter at 
z ~ 0.34 or an [O ill] emitter at z ~ 0. In the case of an 
[On] emitter at z ^ 0.34 we would expect to see [Oiii] 
at A 6710 A, which we do not see. If the object were 
a z [O III] emitter, we would expect to see the A = 
4960.295 A[Oiii] line with a line flux of ~ 1.8 x lO'^'^^ 
which we also don't see. The more likely scenarios are 
then that the [Oiii] emission in this object is weak, or 
the [Oiii] line is being covered by strong OH lines/H2 
absorption in this region. 

We are unable to detect H/3 and/or determine upper 
limits for II/3 emission, likely because the redshift of each 
object has placed the H/3 line under strong OH lines 
and/or under H2O absorption features. We do not de- 
tect the [0 11] (3729.875 A) line by visual inspection in 
either LAE27878 or LAE40844. Determining an upper 
limit for this line by fixing the expected [0 11] wavelength 



based on the redshift measured from the [O ill] line did 
not yield a significant upper limit. 

3.3. Systemic Redshifts and Velocity Offsets Between 
[Oui] and Ly-a 

Using the Gaussian fits described above, we measured 
systemic redshifts from the strongest [O ill] line in the 
two objects with detections, finding redshifts of 3.11170 
± 0.00014 and 3.11879 ± 0.00011 for LAE40844 and 
LAE27878, respectively. 

Measuring the redshift of each object using the Ly- 
a line instead of the [Oiii], yields redshifts of 3.11639 
± 0.00021 and 3.12051 ± 0.00021 for LAE40844 and 
LAE27878, respectively, after corrections for the Earth's 
motion. To correct for the Earth's motion we calculated 
topocentric radial velocitiesH for the two observing lo- 
cations (MMT and LET) for the nights the objects were 
observed at each location. The generally accepted inter- 
pretation of this discrepancy in redshift measurements 
from interstellar emission lines and Ly-a is that there is 
a kinematic offset between the lines caused by a large 
scale outflow, an outflow likely driven by active star for- 
mation. 

Assuming that the emission lines originate from a sin- 
gle redshift, we pin the lines to the redshift of the [O ill] 
line and use this frame to calculate a velocity offset be- 
tween the two lines. We justify using the nebular emis- 
sion to deflne the systemic redshift of the galaxy since the 
[O III] emission originates from H il regions surrounding 
ionizing stars. These regions ought to be at the systemic 
redshift of the galaxy. While the Ly-a initially departs 
from the same regions, resonant scattering, which effects 
Ly-a and not [O ill] , changes the observed location of 
Ly-a emission. 

We derived velocity offset between the 5008.240 A 
[O III] line and Ly-a line based on the central wavelength 
of each line, determined by the best flt asymmetric and 
symmetric Gaussians for the Ly-a and [O ill] lines, re- 
spectively. We find velocity offsets of -1-342 ± 18.3 km 
and +125 ± 17.3 km s"! for LAE40844 and LAE27878, 
respectively. The velocity offsets between the [O ill] and 
Ly-a lines are shown in Figure 2. 

Steidel et al. (2010) note that their redshift determi- 
nations for z ~ 2 - 3 galaxies based on NIR Ha measure- 
ments have an inherent uncertainty of ~ 60 kms"^. This 
estimate is based on repeated observations of the same 
galaxy with their 0.76" slit in different positions. The un- 
certainty is explained as arising from the fact that in each 
measurement they are only measuring the velocity of the 
fraction of the flux that entered the slit. We find that our 
[Oiii] measurements for our z ^ 3.1 galaxies should not 
be subject to such a large uncertainty from this effect due 
to our larger slit width (1") and small galaxy sizes. The 
sizes of our galaxies, from ACS i-band half light radii, are 
1.1 and 1.3 kpc for L AE40844, LAE27878, respectively 
(jMalhotra et al.ir2010[) . The corresponding half-light an- 
gular diameters are still less than half the slit width. The 
1.5" diameter fibers should also minimize flux losses for 
our Ly-a observations, however, we concede that our er- 
ror bars may be lower limits due to such systematics we 
may not be fully taking into account. 

http:/ /fuse. pha.jhu.edu/support/tools/vlsr. html 
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TABLE 1 



Galaxy Characteristics 



LAE14310 



LAE27878 



LAE40844 



u* Magnitude 

Narrowband KPNO MOSAIC[Oiii] Magnitude 
g"*" Magnitude 

'^OIII 

Ly-Q Equivalent Width^-'^ 
Ly-Q Line Flux from Narrowband'^'* 
Upper limit on xray / Ly-a Flux Ratio 
O III line flux (A = 5008.240 A)^ 
O III line flux (A = 4960.295 kf 
Olll velocity widtli(A = 5008.240 A)^ 
O III FWHM (A = 5008.240 kf 
^ofi^set ^ from Lya^ 



25.46 ± 0.29 
22.56 ± 0.11 
24.49 ± 0.13 
3.11043 ± 0.00021 



18.7 



17 
20 
+2.25 
2.51 

3.0 



26.54 ± 0.50 
23.34 ± 0.15 
25.47 ± 0.19 
12051 ± 0.00021 
11879 ± 0.00011 

1+34 
■40 
1.42 
1.63 
6.1 

6.96 ± 0.33 
1.47 ± 0.37 

189.3 ± 10.3 
13.0 ± 0.7 

+ 125 ± 17.3 



1183 
9.41"* 



25.56 ± 0.32 
21.82 ± 0.06 
23.66 ± 0.06 
11639 ± 0.00021 
11170 ± 0.00014 



78 



+8 



36.1 



+2.35 
-2.47 



1.6 

35.48 ± 1.15 
14.82 ± 2.24 
281.1 ± 9.8 
19.3 ± 0.7 
+342 ± 18.3 



^ corrected for Earth's motion 

^ Rest Frame, from Narrowband flux, calculated as (Pjvs 
of KPNO MOSAIC[Oiii] filter and 1265 
g+ band. 



Fg)/(Fg/55A - F]vs/1265A) where 55 A is bandpass 
is bandpass of Subaru g+ filter, Fjvs is fiux in narrowband, Eg is flux in 



Calculated as (i^jvs - 
^ kms~l 



- Fg){c/\1)d\ where Ac central wavelength and dA is bandpass of KPNO MOSAIC[Olll] filter 





Fig. 2. — [Olll] line and Ly-o lines with their corresponding 
best-fit Gaussians. Velocity offset between [Olll] and Ly-o line 
for LAE40844 is 342 km s~^. "Velocity offset between [Olll] and 
Ly-a line for LAE27878 is 125 km s~^. Optical spectra have been 
calibrated using the Ly-o line flux determined from narrowband 
imaging. The feature near + 1400 kms-i inLAE27878 is a residual 
night sky line at ~ 20728.17 / 20729.859 A URousselot et al.|[200Ch . 

3.4. 'Blue bump' in LAE4O8U - Velocity Offset of 
Secondary Ly-a Feature 

LAE40844 has another feature of interest in its optical 
spectrum, namely a smaller, secondary Ly-a peak blue- 
■ward of the systemic velocity of the object. See Figure 3 
for a detailed vie'w of this feature. This feature is fit -with 
an asymmetric Gaussian as described for the main Ly-a 
line in section 3.1, but the constants on the left and right 
sides of the Gaussian are required to be equal (and to be 
equal to the pre-determined constant level also described 
in section 3.1) to ensure that the main Ly-a peak did not 
interfer 'with our best fit measurements of this secondary 
peak. This essentially reduces the number of fitted pa- 
rameters for the blue bump from 5 to 4, meaning that 
"when both the main Ly-a line and the blue bump are 
fit, a total of 9 parameters are returned (5 for main Ly-a 
peak, 4 for blue bump). Our method yields a velocity 
offset from the [Oiii] line of -453.7 ± 50.7 kms"\ after 
correction for the Earth's motion. From this measure- 
ment "we determine that the t"wo Ly-a peaks are offset 
from one another by +796.2 ± 53.9 km s~^. 




-1500 -100c -500 500 1000 '500 



Fig. 3. — LAE40844, best fit asymmetric Gaussian to redshifted 
Ly-a emission in red, best fit asymmetric Gaussian to blue shifted 
Ly-a emission in blue, optical spectrum in black. Velocity offset 
between the Ly-a two peaks is -p796.2 ± 53.9 km s~^. 

Additionally, using the flux calibration "we derived from 
the narro'wband line flux, we find that this blue Ly-a 
peak has a line fiux of ^ 1.08 xlO^^^ erg s~^ cm^^. 
When determining this calibration "we found that based 
on the transmission curve of the narro"wband filter, the 
blue bump contributed ~ 9.4 % of the total narro'wband 
line flux (3.61 xlO~^^ erg s~^ cm~^). Comparing this 
to a line flux of ~ 3.27 xlO~^^ erg s~^ cm~^ for the red 
Ly-a peak "we find an approximate flux ratio, red:blue, 
for the t'wo lines of 3.0. In other 'words, the strength of 
secondary (blue) peak is roughly 33% that of the main 
(red) peak. In section [5] 'we discuss a scenario that can 
give rise to this blue bump and compare the velocity off- 
set we find bet'ween the t'wo Ly-a lines to velocity offsets 
that have been presented in the literature on Ly-a radia- 
tive transfer. 

3.5. Effect of the [Olll] Emission Line on Mass and 

Estimates 



The COSMOS field has a deep K^-band (centered at 
21460A, AA = 3250 A) coverage fromCFHT WIRCAM, 
'with a 5(T depth of ^ 23.8 mag in a 3" aperture on a 
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PSF matched image. For our z ~ 3.1 LAEs, this band 
win encompass both the continuum and [O ill] emission. 
The 2008 COSM OS Intermediate a nd Broad Band Pho- 
tometry Catalog (jCapak et al.ll200'7l ) has Kg band magni- 
tudes of 22.61 ± 0.07 for LAE40844, and 24.91 ± 0.62 for 
LAE27878 (MAG_AUTO measurements from a 3" aper- 
ture) . Our measured [O ill] hne fluxes can account for the 
entire Kg band fluxes, where we find K^-band magnitudes 
of 23.09 ± 0.036 and 24.79 ± 0.053 for LAE40844 and 
LAE27878, respectively, using just the [Oiii] line fluxes. 
This shows that our detected [Olll] emission lines alone 
can be responsible for all or nearly all the flux measured 
in the Kg band for both of these LAEs. 

This has important implications for mass and age es- 
timates of high redshift galaxies. These estimates typi- 
cally rely on the size of the Balmer break to determine 
the age of the galaxy. If, as is the case for our two 
LAEs, there is a dominant emission line polluting the 
location redward of the Balmer / 4000 A break, then the 
size of the break may be overestimated and the subse- 
quently derived ages and stella. r mass es may be overes- 
timated. iSchaerer k, de BarrosI (J2009) found that when 
nebular em i ssion l ines were included, SED fitting of the 
lEvles et all ()2007[ ) sample of 10 z '--^ 6 galaxies yielded 
an average age ~ 4 times younger than what was found 
without the emission lines included. The average stellar 
mass estimate also decreased from 1.2x10"'^° M© with- 
out emission lines to 7.8x10^ M© when nebular emission 
was included. Some studies of Ly-a emission in Lyman- 
break selected populations have found that the strongest 
Ly-a emitt ers have blue UV sp e ctral slopes but red opti- 
cal slopes (jShaplev et al.ll2003l: iKornei et al.ll2010t ). and 
have concluded that Ly-a emission is strongest in LBGs 
that are older but relatively dust-free. If strong nebular 
line emission contributes to the observations of red rest- 
optical slope, it might be possible to reinterpret such 
observations in te rms of young , stron gly line emitting 
galaxies, although IKornei et al.l (|2010D argue that their 
observed correlations between Ly-a strength and stel- 
lar population age are unchanged when they select ob- 
jects only at redshifts where the optical continuum filters 
are line-free. Previous work on stellar populations of 
Ly-a-selected galaxies has found that nebular line emis- 
sion is required to explain observed rest-optical colors 
(jFinkelstein et al . 2008). Recently, at the highest red- 
shifts. [Ono et al i ((2010) have shown that either old stel- 
lar populations or young ones with strong nebular emis- 
sion can reproduce the composite SEDs of Ly-a selected 
galaxies. The older models have correspondingly higher 
stellar masses, since mass to light ratio increases strongly 
with age. Our observations provide direct observational 
evidence that nebular line flux dominates the rest optical 
in analogous objects at z sa 3, and hence supports the 
interpretation of the high redshift Ly-a selected popula- 
tions as young and low-mass objects. 

4. DISCUSSION 
4.1. [Olll] Luminosities 

The [O III] line has been measured in other objects 
at similar redshifts. Pettini et al. (2001, henceforth 
POl) find [Oiii] (5007 A) luminosities of 3.5 - 15.6 x 
-|^q42 gj.g g-1 individual z ^ 3.1 LBGs observed 

with VLTl/ISAAC and/or Keck II/NIRSPEC. Compar- 



ing this to our range of luminosities for the 5007 A line 
for two z - 3.1 LAEs of 6.1 - 31.0 x lO'*^ erg s-\ it 
appears that our fainter 5007 A measurement falls in the 
POl range, while the stronger of our two 5007 A lines 
is almost double that of the brightest luminosity in the 
POl sample. Luminosities for the 4959 A line in the POl 
sample are 1.4 - 6.5 x 10''^ erg s~^. This yields the same 
trend we see in the 5007 A line; where our fainter 4959 
A measurement falls in the POl range, and our brighter 
4959 A line is approximately twice that of the brightest 
4959 A measurement in the POl sample. 

[O III] luminosities from Ic nsed galaxies a r ound z ^ 3 
have also been documented. iFosburv et al.l (|2003[ ) mea- 
sured the [O III] line in a lensed Hl i galaxy, also known 
as the Lynx arc (jHolden et al.l[200ll ) , at a redshift of z ~ 
3.36. Using NIRSPEC K-band spectra they find an [O ill] 
luminosity of 28.3 ± 0.3 x 10^^ gj.g for the 5007 A 
line and 9.8 ± 0.3 x lO'*^ erg s'^ for the 4959 A hne. 
These luminosities (which have been corrected for mag- 
nification) are quite comparable to our measurements for 
the 5007 A and 4959 A lines. iFinkelstein et al.l (|2009bf ) 
measured the [O ill] line in a lensed ultraviolet-luminous 
z = 2.73 galaxy known as the 8 o'clock arc. For the 4959 
A line they find a luminosity of 8.9 ± 0.4 x 10'*^ erg 
(after correction for magnification) . This again falls right 
in the range reported for our two LAEs . The 50 07 A line 
was not reported by IFinkelstein et al.l (|2009b( ) because 
it fell in an area of low atmospheric transmission which 
required a correspondingly large tellluric correction and 
led to large uncertainties in any measurements from the 
line. 

Looking at our sample of two galaxies, the [O ill] lumi- 
nosities in LAE27878 are most similar to the POl LBGs 
while the more luminous [O ill] lines in LAE40844 are 
more comparable to more luminous [O ill] lines found in 
galaxies studied via lensing. A larger sample of NIR spec- 
tra with [O III] line measurements for LAEs will help us 
understand the range of [O ill] luminosity in LAEs and 
its implications. 

4.2. Ly-a, Line Profiles and Outflow Models 

In addition to the information that can be gleaned from 
the line fluxes of nebular emission lines such as [O ill] , 
the asymmetric profiles of Ly-a emission lines themselves 
carry information on the physical conditions and pro- 
cesses in these objects. Resonant scattering can lead to 
asymmetric profiles through radiative transfer processes 
operating either within the Ly-a emitting galaxy, or in 
the surrounding intergalactic medium. Establishing the 
systemic velocity with the [O ill] line sheds new light on 
those processes and conditions. 

We find that the Ly-a line profiles seen in LAE40844 
and LAE27878, and the velocity offset of the Ly-a line 
from systemic in both objects is in good agreement with 
what is predicted by outflow models, where the Ly-a 
is redshifted through interaction with receding gas on 
the far side of the galaxy, and transmitted through ap- 
proaching gas on the near side due to the line's kine- 
matic redshift. Two particular types of outflow models 
are discussed here due to their apparent agreement with 
our resul ts - the ca se where the outflow is in a coher- 
ent shell (jTenorio-Tagle et al...l999; .Dawson et al.l 120021: 



iZheng et aLll2010l: iVerhamme et all 
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2006 



1991 



l2008l) and the 
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ansen 



fcOhl 



case of a clumpy outflo w (jNeufeL 
[2001: ISteidel et alJ[20lol) . 

Understanding our results in the context of an out- 
flow is justified in that high-redshift LAEs are typically 
youn g (age ~ 10^ years) wi t h vigorous star formation 
(e.g.. iFinkelstein et all [20071 I2Q09I: iPirzkal et alj 120071: 
iGawiser et al.ll2007t iLai et al.ll2008l ). Their typical star 
formation intensities are well a bove the threshold re - 
quired to drive galactic winds (jMalhotra et al.' '201(f). 
Similar winds are seen in nearby starburstin g galaxies, 
with velocities of ord er 10^-10^ km s^^ (Hcc kman et al.l 
ll990HHeckniaiill2002D . numbers that encompass our mea- 
sured velocity offsets. 

In the model with a single expanding shell, with a 
central monochromatic source (Verhamme et al. 2006, 
henceforth V06), the redshifted Ly-a line is built up of 
photons that underwent one or more backscatterings off 
the expanding shell. The more backscatterings the pho- 
ton undergoes, the further it is redshifted, giving rise to 
the prominent red wing that is seen in the redshifted Ly- 
a line. Photons that are emitted from the blue wing of 
the Ly-a line in the part of the shell that is approaching 
the observer can give rise to the blue bump we see in 
LAE40844. See Figure 12 of V06 for a detailed descrip- 
tion of photons escaping from an expanding shell. 

For the parameter space examined in the V06 simula- 
tions, velocity offsets of a few 100 kms^^ are predicted 
for the redshifted Ly-a line. The magnitude of the Ly-a 
velocity offset depends on the velocity of the expanding 
shell, the column density of neutral hydrogen and the 
Doppler parameter (see V06 for details on these parame- 
ters). Our velocity offsets of 125 kms~^ and 342 kms~^ 
appear to be at the low end of this range. When the 
blue bump is seen in these simulations it is offset from 
the redshifted Ly-a line by ^ 1000 kms^^, in agreement 
with our measured offset in LAE40844 between the red 
and blue peaks of 796 kms^^. Thus there is good gen- 
eral agreement between our observations and the V06 
models, which have not been tuned specifically to fit our 
data. It seems likely that an expanding shell model could 
fit our data well with some adjustment of the input pa- 
rameters. Additionally, while our total offset between 
the blue and red Ly-a peaks in LAE40844 agree with 
the total offsets seen in their simulations, we find that 
^biuepeak ~ -2vredpeafc in LAE40844. lu V06 the velocity 
offset between these two peaks is nearly symmetric in the 
cases where both peaks are present and the expansion ve- 
locity of the shell is small (< 200 kms~^). Where the 
expansion velocity of the shell is large (300 - 400 km s^^) 
the blue peak is nearly -1/2 the velocity shift of the red 
peak. In either case, these predictions do not directly 
match our observations presented here. Further work is 
needed to understand the discrepancy between the ve- 
locities of the redshifted and blueshifted peaks in models 
and observations. Deviations from spherical symmetry 
of the expanding shell mo del could help account for this 
difference, as was noted in lSchaerer fc Verhammei (j2008[ ). 

These same authors further investigate their model in 
a later paper (Verhamme et al. 2008, henceforth V08) by 
applying it to fit the observed Ly-a emission line profiles 
of the LBGs from Tapken et al. (2007), which lie at 
2.7 < z < 5. Varying model parameters include the 



Doppler parameter, Yexp, the neutral hydrogen column 
density Nhi , the dust optical depth Tq , and the intrinisc 
EW and FWHM of the input Ly-a line. Nine of the 11 
LBGs investigated were found to have Yexp ~ 150 - 200 
kms~^. Two LBGs, which resemble our LAE40844 in 
that they have a blue bump, can be fit either with larger 
velocities, of order 300 - 400 kms~^ or a with a quasi- 
static medium with Vexp of order 10 - 25 kms""'^. V08 
prefer the quasi-static explanation. We note that our 
results oivoffset — 125 and 342 kms^^ is fully consistent 
with the V08 results. However, we emphasize that we 
measure the velocity offset between the Ly-a line and 
the systemic redshift as defined by [Oiii], whereas V06 
and V08 are quantifying Yexp, which is the velocity of 
a spherically expanding shell around the central Ly-a 
source in their model. We also note that LAE40844, 
with its blue bump and Voffset of 342 kms~^, appears 
to show that a quasi-static model is not always preferred 
for objects with a blue bump. 

Steidel et al. (2010, henceforth SIO) consider out- 
fiows to be an important component of the mechanism 
that shapes the observed Ly-a profile, but they prefer 
a scenario in which the structure and kinematics of the 
circum-galactic medium can produce our observed pro- 
files, instead of relying on Ly-a radiative transfer in an 
expanding shell to create the profiles we see. In the SIO 
scenario, a clumpy outfiowing gas will allow some Ly-a 
photons to escape from a galaxy producing the redshifted 
Ly-a line we have observed. This same scenario is also 
capable of producing the faint 'blue bump' we have dis- 
cussed in section 3.4. SIO find an average velocity offset 
of 445 kms~^ in a sample of 42 z ~ 2 - 3 LBGs with 
Ly-a emission. The velocity offset measurement in this 
case was made with respect to the redshift defined by 
Ha emission. As is the case when we compare our re- 
sults to the Verhamme model, our velocity offsets of 125 
kms""'^ and 342 kms~^ for our observed LAEs are below 
these values. This conclusion includes consideration of 
the fact that SIO measures the centroid of a single func- 
tion whereas we measure the peak of Ly-alpha through 
an asymmetric function composed of two Gauss func- 
tions; which inherently causes the SIO measurements to 
be higher. This causes a 10 - 15% difference in the veloc- 
ity offsets calculated, and hence we stil find our velocity 
offset values to be lower than those in the SIO sample 
evn after this consideration. 

At present, our sample of two z ^ 3.1 LAEs reported 
in this paper does not allow us to distinguish between 
the expanding shell scenario or the SIO interpretation, as 
both are able to produce profiles and velocity offsets in 
reasonable agreement with our observations. Whichever 
of the scenarios discussed here (V06 or SIO) is pro- 
ducing the observed velocity offsets, we emphasize that 
winds/outflows are important in either case. 

The additional observational samples to which we can 
best compare our current results are z 3 LBGs. POl 
and Shapley et al. (2003, henceforth S03) have both mea- 
sured the velocity offset of the Ly-a line from systemic 
in a population of z 3 LBGs. POl find velocity offsets 
for the Ly-a line of 200 - 1100 kms~^ in a sample of 
13 LBGs that also show Ly-a in emission. The velocity 
offset is measured compared to the redshift of nebular 
Hii emission. S03 find a velocity offset of 360 kms~^ 
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from a composite spectrum of 811 z 3 LBGs. The 
Ly-a velocity offset measurement in 803 was made with. 
respe ct to interstehar absorption hnes. (jTapken et aLl 
l2007f l measured a velocity offset between LIS lines and 
Ly-a in seven LBGs with Ly-a emission at redshifts of 2.7 
- 5 (the same sample of galaxies used in V08) finding an 
average offset of 580 kms^^. Including the 445 kms~^ 
offset measurement from the SIO sample discussed above, 
we note that in all cases our observed velocity offsets of 
125 - 342 kms^^ in two z ~ 3.1 LAEs arc at the low 
end of the values reported for the various LEG samples, 
although the velocity offset measurements have thus far 
been made via different methods. 

An alternative explanation of the observed Ly-a pro- 
files, not based on galactic scale outflows, is that they 
arise through resonant scattering in the intergalactic gas 
surrounding an LAE. Zheng et al. (2010) have explored 
such a mechanism in detail for redshift z = 5.7. Their 
models can produce Ly-a lines that qualitatively resem- 
ble our observations both in the line asymmetry and in 
the redshift of the Ly-a line. However, some caution is 
needed in applying these results to our data set, given 
that the IGM density at redshift z = 3.14 is ^ 1/4 that 
Sit z — 5.7, and the ratio of neutral gas density between 
these two redshifts is still more extreme. Overall, we con- 
sider it more likely that winds play an important role in 
Ly-a escape, given that winds are generically expected 
for galaxies with the high specific star formation rates 
typical of LAEs. 

4.3. Implications of Detected Outflows 

Detection and characterization of galaxy scale out- 
flows at high redshift is important because these out- 
flows have important consequences for the evolution 
of individual galaxies as well as the evolution of the 
IGM. Large scale galactic outflows are capable of driv- 
ing materials out of the galaxy and may therefore con- 
tribute to metal enrichment of the IGM at high-z by 
introducing materials produced from starbursts into the 
IGM. In addition, galactic winds likely provide a cru- 
cial channel b y which ionizing photons can escape from 
a galaxy (Steidel et aI]|2001l:IHeckman et al.|[2b01,) . This 
has important implications for the contribution of high- 
z galaxies to the reionization of the IGM. In terms 
of shaping an individual galaxy, superwinds are re- 
sponsible for driving d ust from a starbursting galaxy 
(iFerrara et al.lll991: Heckmaii et al.| [2000: Shaplc v et al.l 
120011: IBianchi fc Ferrarall2005ir and the mass loss from a 
galaxy due to a n outflow may be capable of suppressing 
star formation (jSomerville fc Primackl 119991 : iHeckmanI 
[200^ . While our work presented in this paper has now 
demonstrated that that our sample of two Ly-a selected 
galaxies at z ~ 3.1 are driving winds, further character- 
ization of these winds from a larger sample will help us 
understand and test some of the broader implications of 
winds detailed above. 

The observed velocity offsets between the Ly-a line 
and systemic velocity also have important implications 
for L y-a based t es ts of r eio nization ( Malhotra &: Rhoads 
200l ISantosI ^OOl IMalhotra fc Rhoadsl 12001 
Diikstra fc Wvithd I2010D . In particular, the ion- 
ized volume test proposed by IMalhotra fc Rhoadsl 
(|2006f) (hereafter MR06) is much more sensitive if Av 
(Ly-a peak vs. systemic) is typically small. That test 



works by noting that substantial transmission of Ly-a 
through a generally neutral IGM requires a locally 
ionized region around each observed LAE. The product 
of the bubble volume V and the LAE number density 
n is then a filling factor of ionized gas, from which the 
volume fraction of the ionized phase is « 1 — exp{—nV) 
(MR06). The relevant bubble volume V is sensitive 
to the velocity offset, since Ly-a photons that are 
already redshifted before leaving the emitting galaxy 
are less strongly scattered by the damping wing of 
neutral hydrogen in the surrounding intergalactic gas. 
The effect is explored in detail in Figure 1 of MR06. 
Replacing the range < Av < 360 kms~^from that 
paper with our measurements, Av — 125 kms~"'^and 
At; = 342 kms""'^, would narrow the range of permitted 
volume ionized fractions from the 20% - 50% range 
derived in MR06 to - 20% - 40%. While this discussion 
is subject to refinement as the sample of LAEs with a 
measured velocity offset grows, it shows the importance 
of measuring Av for studying reionization with Ly-a 
lines. 

Finally, even at redshifts where the IGM is predom- 
inantly ionized and affects the Ly-a line only through 
the Ly-a forest, our systemic redshift measurements 
have important implications. Several groups have shown 
that the observed Ly-a luminosity func tion is largely 
unchanged from z « 3 to 2: 6 (e.g., 'Dawson et ajj 

t34; Ster n et al.l l2005: Oui chi et al.ll2 008: Cassata ct al| 
11). Recently, Cass ata et al.l (12011 ) have combined 
this observation with the e xpec ted o ptical dept h evolu- 
tion of the Ly-a forest (Madaul ll99"5l : iFan et aLll2008D to 
infer that the Ly-a luminosity function is in fact evolv- 
ing towards higher luminosities at higher redshifts. The 
key assumption in this argument is that the line emitted 
by the LAE is symmetric and centered on the systemic 
velocity, so that the fraction scattered by the IGM ap- 
proaches 50% by z « 6. This implies that the fraction 
of Ly-a fiux observed would decline by a factor of « 0.6 
between redshift z « 3 and z ~ 6. However, for the 
two objects where we observe [Olii], we know that only 
those photons observed blueward of the systemic velocity 
would be subject to additional Ly-a forest absorption at 
higher redshift. For LAE40844, the blue bump would be 
progressively obscured by the forest at higher redshift, 
resulting in a flux loss of a factor « 3.0/4.0 = 0.75 or 
so at most. For LAE27878, there is no significant flux 
blueward of the systemic velocity, and the fraction of 
its Ly-a emission that we can see should remain nearly 
unchanged from z = 3.1 until the IGM neutral fraction 
becomes so large that the red damping wing of the IGM 
becomes optically thick — i.e., until we reach the central 
stages of reionization. 

To recap, our objects show relatively little Ly-a flux 
blueward of the systemic velocity. This either requires 
dropping the assumption that Ly-a is symmetric and 
centered at the systematic velocity, or else implies that 
we happen to have observed two sources where the Ly- 
a forest at z 3.1 is unusually optically thick. We now 
proceed to evaluate the likelihood of the second scenario. 

Presently, our sample consists of two galaxies, with 
measured blue to red flux ratios of 0.33 and 0.08 for 
objects 40844 and 27878, respectively. We estimate 
the chance of such an occurrence under the assump- 
tion of Cassata et al (i.e., that the intrinsic ratio is 1:1, 
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and deviations indicat e absorption by the Ly-a forest). 
iMcDonald et al.l ()2000D give probability distributions for 
transmission through the Lyman a forest at z = 3.00 
and z = 3.89. Interpolating their results to our redshift 
(z = 3.12), we infer probabilities P{T < 0.33) = 0.24 
and P{T < 0.08) = 0.15. Given a sample of 2, and treat- 
ing their Ly-a forest transmissions as independent ran- 
dom variables, the probability that at least one will have 
T < 0.08 while the second has T < 0.33 becomes 0.050. 
So, our present results disfavor this assumption of an in- 
trinsic 1:1 ratio of blue:red flux, suggesting t hat th e lumi - 
nosity function evolution inferred in Cass ata et al.l (|2011[ ) 
is a consequence of their implicit assumption Av = and 
not a strong conclusion about the true evolution of Ly-a 
galaxy populations. For now, this is a 2a result. Ob- 
servations of a few more Ly-a emitters with systemic 
velocity measurements could resolve this question firmly. 



5. CONCLUSIONS 

We have detected [O III] emission in two Ly-a selected 
galaxies at z ~ 3.1 using the new NIR spectrograph LU- 
CIFER on the LBT. This is a successful demonstration 
that the [O ill] line can be detected in high-z Ly-a se- 
lected galaxies and that this line can be used to investi- 
gate the characteristics of these galaxies. 

In both LAEs we measured a velocity offset between 
the Ly-a emission and the systemic redshift of the galaxy 
as defined by the [O iil] emission. These velocity offsets 
range from 125 - 342 kms~^. We find that these veloc- 
ity offsets and the observed profile of the Ly-a line both 
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